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MARKED-UP TRANSLATION OF INTERNATIONAL APPLICATION 

D e scr i pt i on TITLE OF THE INVENTION 

METHOD AND ARRANGEMENT FOR DESIGNING A TECHNICAL SYSTEM 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to PCT Application No. 
PCT/DE2003/002566 filed on July 30, 2003 and German Application No. 102 37 335.3 filed on 
August 14. 2002, the contents of which are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

[0002] The invention relates to a method and an arrangement for designing a technical system. 

[0003] In order to design a complex technical system it is often necessary to optimize saidthe 
system i n r e sp e ct o f with respect to a plurality of contradictory criteria. Sa&The criteria influence 
target functions of the system, such as, for example, manufacturing costs or efficiency. In 
addition, possible operating points of the system can be restricted by auxiliary conditions. This 
leads to the problem of determining a set of optimal operating points for the system, that is to 
say the set of possible operating points of the system with which it is not possible to optimize 
saidthe operating points further simultaneously with regard to all criteria. From saidthe set of 
optimal points A individual users can then select the most suitable operating points of the system 
for their applications while taking into account secret criteria or expert knowledge. 

[0004] A weighting method for optimizing technical systems mwith respect ofto a plurality of 
criteria is known from ffl C. Hillermeier: "Nonlinear Multiobiective Optimization: A Generalized 
Homotopy Approach", Chapter 3.2. Birkhauser Verlag. 2001 ("the Hillermeier Chapter 3.2 
reference") , wherein scaling parameters are employed to apply transformations to scalar-value 
optimization problems. Sak lThis method has the disadvantage that it is numerically very 
involved, because very many scalar-value optimizations have to be performed. Furthermore, the 
selection and variation of the scaling parameters necessitates an interaction with a user and in 
this respect cannot be automated. 

[0005] A stochastic method for optimizing technical systems mwith respect ofto a plurality of 
criteria, wherein a stochastic differential equation is used to solve the optimization problem, is 
described in fSi C. Hillermeier: "Nonlinear Multiobiective Optimization: A Generalized Homotopy 
Approach", Chapter 3.3. Birkhauser Verlag. 2001 . Sate This method has the disadvantage that it 
i eare very involved in numerical terms, because a multiplicity of quadratic optimization problems 
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have to be solved. A further disadvantage lies in the fact that with ea&the method A the individual 
target functions are not weighted, as a result of which important information for selecting an 
optimal point is not available to the user. 

[0006] A homotopy method for optimizing technical systems m with respect efto multiple criteria, 
wherein in addition to weighting factors for the target functions A Lagrange multipliers are used in 
order to take auxiliary conditions into account, is known from \3\ C. Hillermeier: "A Generalized 
Homotopy Approach to Multiobiective Optimization", Journal of Optimization Theory and 
Application, Vol. 110/3, pp. 557-583. Plenum Press, New York. 2001 ("the Hillemermeier Vol. 
110/3 reference") . The disadvantage of saidthis method lies in the fact that an interaction with 
the user is necessary and in this respect satethe method cannot be automated. 

SUMMARY OF THE INVENTION 

[0007] Tfr eOne possible object of the invention is therefore to create an automated and 
numerically efficient method for designing a technical system. Th i s obj e ct is ach ie v e d according 
to th e f e atures of th e i nd e p e nd e nt c l a i m s . D e v el opm e nt s of th e i nv e ntion w i l l also b e appar e nt 
from th e d e p e nd e nt claims. 

[0008] W i th th e The inventors propose a method accord i ng to th e i nv e nt i on f or designing a 
technical system in which 7 the technical system is modeled by m e ans of by a predetermined set 
of target functions which are dependent on parameters. In this modeling process each individual 
target function is weighted with a weighting factor. The method solves an equation system 
comprising the parameters and the weighting factors as variables in a variable space, with 
solutions of the equation system forming operating points of a solution space in the variable 
space. In the method accord i ng to th e inv e nt i on the operating points are determined by m e ans 
efby a predictor-corrector method, according to which, starting from a first operating point, a 
predictor generated by m e ans of by a stochastic variable is determined in the variable space, 
and subsequently, in a corrector step, a second operating point is determined. SaidThe 
determined operating points are used here to design the technical system. The i nv e nt i v e 
method can be used to design a new technical system or modify or, as the case may be, adapt 
an existing technical system. 

[0009] An advantage of th e i nvention resides in the fact that the method is automated through 
the generation of the predictor by m e ans of by a stochastic variable and so there is no longer 
any need for afv-intervention on the part of the user. The linking of the numerical predictor- 
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corrector method with stochastic methods guarantees efficient use of the computer resources 
for calculating operating points of a technical system. 

[0010] In an advantageous embodiment of th e i nv e nt i on the predictor is predetermined by 
m e ans of by random numbers, so that in particular during the execution of the method a random 
number generator can be used and b vthrough this A means-the automation of the method is 
ensured in a simple manner. 

[0011] In a further particularly advantageous embodiment of th e inv e ntion the random numbers 
are normally distributed. What this achieves is that the trajectory of operating points which forms 
in the solution space during the execution of the method ensures a uniform distribution in the 
entire solution space and so ensures good coverage of all possible operating points. As a result 
of the use of normally distributed random numbers, in particular a Brown mov e m e nt Brownian 
movement on the solution space can be modeled by m e an6 of by the method accord i ng to th e 
inv e nt i on . 

[0012] Preferably the operating points which are determined by m e an s of by the method 
accord i ng to th e i nv e ntion are what are known as pareto-optimal points which cannot be 
optimized further in relation to all target functions. In the method accord i ng to th e i nv e nt i on , in 
particular the points with positive weighting factors in the solution space are selected as 
operating points. 

[0013] In a further advantageous embodiment of the inv e nt i on the operating points must also 
satisfy one or more auxiliary conditions, with the or each auxiliary condition being represented 
by a further variable of the equation system in the variable space. In this case the auxiliary 
conditions can be equality auxiliary conditions and/or inequality auxiliary conditions. With 
inequality auxiliary conditions a slack variable is preferably introduced, by m e ans of by which the 
inequality auxiliary conditions can be transformed into equality auxiliary conditions. The use of 
slack variables will be explained in more detail in the detailed description of an exemplary 
embodiment. 

[0014] The solution space of the operating points is preferably a manifold, in particular a 
submanifold in the variable space. In te tthe Hillemermeier Vol. 110/3 reference it is explained 
under what preconditions the solution space forms such a manifold. 
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[0015] Since in particular at the start of the method according to tho i nv e ntion a first operating 
point is present initially, in a special embodiment of th e i nv e ntion t his first valid operating point is 
determined by m e ans ofby a weighting method, the use of weighting methods already being 
known from the pfioj^af trelated art (see fflthe Hillermeier Chapter 3.2 reference) . 

[0016] With the predictor-corrector method , which may be used in the invention, a tangential 
plane to the solution space is determined, preferably in the first operating point, and the 
predictor is then specified in saidthe tangential plane. 

[0017] In a development of the method accord i ng to th e i nv e ntion , if a negative predictor with 
one or more negative weighting factors occurs, a new predictor is determined by a reflection at 
a subplane of the solution space of the valid operating points. By- Through t his, means-new 
regions of valid operating points can be determined, which operating points can be of particular 
relevance to the user in terms of secret supplementary criteria or his/her expert knowledge. 

[0018] In a preferred embodiment of th e i nv e nt i on , in the reflection step a point of intersection of 
the trajectory that runs between the first operating point and the negative predictor with a 
subplane of the solution space is determined. The tangential component of the vector spanned 
by the point of intersection and the negative predictor to the relevant subplane of the solution 
space is then determined, with those weighting factors which were negative for the negative 
predictor in the points of the subplane now being equal to zero. Next, the normal component, 
associated with sai4the tangential component, of the vector spanned by the point of intersection 
and the negative predictor is determined. Finally, the new predictor is determined by two times 
subtraction of the normal component from the negative predictor. 

[0019] A Newton method known from the pfio^af trelated art , which method is easily convertible 
numerically, is preferably used for the corrector method e mp l oy e d in th e i nvent i on . 

[0020] The operating points are preferably determined by m e ans of by iterations of the predictor- 
corrector method, with the second operating point of the preceding iteration step being used in a 
current iteration step as the first operating point of the predictor-corrector method. In this case 
the method is terminated by m e ans of by, for example, an abort condition. In an advantageous 
embodiment of th e i nv e ntion the abort condition is met when a predetermined number of 
operating points has been determined and/or a predetermined time limit has been reached. 
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[0021] In addition to the above-described method for designing a technical system, the 
inventors propose i nv e nt i on a l6 o r el at es to an arrangement for designing a technical system by 
m e ans of by which the above-described method can be performed. In particular the method 
comprises a processor unit by m e ans of by which it is made possible for the predictor to be 
generated using a stochastic variable. 

[0022] The arrangement preferably comprises a random number generator for generating 
random numbers which represent the stochastic variable. 

[0023] Th e i nvention or any d e v el opm e nt th e r e of d e scr i bed i n th e for e go i ng can also b e 
i mp le m e nt e d by m e an 6 of The inventors also propose a computer program product which has a 
storage medium on which is stored a computer program which is executable on a computer and 
executes the inv e nt i on design method . 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0024] Ex e mp l ary e mbod i m e nts of th e inv e nt i on ar e d e scr i b e d and exp l a i n e d i n th e fo ll owing 
with r e fer e nc e to th e drawings, i n wh i ch: These and other objects and advantages of the present 
invention will become more apparent and more readily appreciated from the following 
description of the preferred embodiments, taken in conjunction with the accompanying drawings 
of which: 

Fig. 1 shows a flowchart of the method according to one embodiment of t he invention for 
designing a technical system; 

Fig. 2 shows a diagram which illustrates the predictor-corrector method used in one 
embodiment of the invention; 

Fig. 3 shows a diagram which illustrates the reflection method used in a var i at i on an 
alternative embodiment of the invention, and 

Fig. 4 shows a processor unit for performing the method according to one embodiment 
ofthe invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0025] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference numerals refer to like elements throughout. 
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[0026] A flowchart of toea method accord i ng to th e i nv e ntion for designing a technical system is 
shown in Fig. 1. 

[0027] First, in step 101, a description form of the technical system is selected. Sa&The 
description form comprises a predetermined number k of target functions f=(f1,...,fk), with each 
of the target functions being dependent on n predetermined parameters x-| to x n of the technical 

system. The target functions are for example the investment costs f<| and the efficiency of the 

technical system. In this case the target functions are described by m e ans of by the following 
equation: 

r , v , Investment. cos ts(x) s 

J (*) - ( — —) = ( — ) 

~ f 2 (x) _Efficiency(x) 

where x=(x p ..., x„). 

[0028] The parameters x-| to x n can be configuration parameters or operating parameters of the 
technical system. 

[0029] By m e an s of By the method accord i ng to th e i nv e nt i on t he valid operating points which 
are used for the design of the technical system are determined by the optimization of the target 
functions i n r e sp e ct o f with respect to the parameters, whereby not all target functions f 1 to f^ 

can be optimized simultaneously since the optimization criteria are generally in competition with 
one another. 

[0030] The technical system is further limited in the valid operating points by a predetermined 
number m of auxiliary conditions h=(hi(x)...,h m (x)) which can be expressed by m e ans of by the 

following equation: 
h(x)=Q 

where 0=(0,...,0) represents a zero vector. In this case what is involved is an equality 
auxiliary condition, with an inequality auxiliary condition also being able to be considered as an 
alternative. An inequality condition of this kind is, for example: 

h(x)<Q or h(x)>Q . 
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[0031] In order to solve the optimization problems by m e an6 of by inequality auxiliary conditions, 
a number m of slack variables s=(si,...,s m ) are introduced, by means ofby which the inequality 

auxiliary conditions can be transformed into the following equality auxiliary conditions: 
h(x)+s=0 or h(x)-§=0 

[0032] With the optimization method described in the present embodiment of the i nv e nt i on , the 
valid operating points are known as pareto-optimal points, which satisfy the following condition: 

min^ {fW},R = {xeM n \h(x) = 0} 

[0033] It can be shown that the solutions of this optimization problem are the solutions of the 
following nonlinear equation systems: 



'XL a.-V/ 1 (x) + Z7 -i A / V* y (x) > | 
h(x) 



= 0 



[0034] In this case the auxiliary conditions are taken into account by the Lagrange multipliers 
hr&V—&m*l and the tar 9 et functions fj are weighted with weighting factors aj, whereby care 
must be taken to ensure that the total of all weighting factors is normalized to one, 
i.e. Zt/ a, - 1 = 0 . In this case, in particular, the individual weighting factors aj can also be 

negative or equal to zero. The solutions of the optimization problem are therefore vectors (x,^,a) 
in the (n+m+k)-dimensional variable space of the above equation system. 

[0035] As shown in43 1 the Hillemermeier Vol. 110/3 reference , under certain conditions the 
solutions of this equation system describe a (k-l)-dimensional submanifold M in the variable 
space. 

[0036] The below described numerical steps for determining valid operating points are 
essentially based on the homotopy method described in43i the Hillemermeier Vol. 110/3 
reference , wherein a predictor-corrector method is used for determining pareto-optimal points. 

[0037] In step 102, proceeding from the description form 101 of the technical system, a first 
pareto-optimal point z is determined by m e an s of by a standard method such as, for example, 
the weighting method. 
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[0038] In this first pareto-optimal point, in the next step 103, a (k-l)-dimensional tangential 
plane T 2 M to the manifold M of the valid operating points is determined in point z. Toward that 

end, a Jacobi matrix of the equation system F in point z is subjected to a QR factorizing. From 
this, an orthonormal basis { q1...q[<_i } is then determined which spans the tangential plane. The 

individual numerical steps performed in this process are described in detail in431 the 
Hillemermeier Vol. 110/3 reference . 

[0039] In the next step 104, a predictor y is determined in this tangential plane, with the 
predictor - in contrast to the homotopy method described in f3i the Hillemermeier Vol. 110/3 
reference - being generated by m e an 6 ofby a normally distributed random number vector b of 
the dimension k-1 in the tangential plane. In this case the predictor y has the following form: 

y =z +(qi...q k _i)b 

[0040] Through the use of a random number vector such as this, a Brown mov e m e nt Brownian 
movement can be modeled on the submanifold M, with the Brown mov e m e nt Brownian 
movement being able to be represented approximately as follows: 

dZ t =eP(Z t )dB t 
where 

P(z) is a projection matrix onto the tangential plane T Z M in the valid operating point z, 
8 is a scaling factor, and 

B t , t € 5Rq is a Brown mov e m e nt Brownian movement in the variable space. 

[0041] In order to model this movement, the k-1 -dimensional normal distribution N(0|<_1 , tAel|<_i) 
is selected for b, where the mean value 0|<_i is the (k-l)-dimensional zero vector and the 
variance is the (k-l)-dimensional identity matrix l|<_<j multiplied by a step increment t^ of the 
Brown mov e m e nt Brownian movement and the scaling factor e. 

[0042] An alternative method of determining the predictor is first to determine a normally 
distributed random number vector in the (m+n+k)-dimensional variable space and then to 
project saidthe vector into the (k-l)-dimensional tangential plane T Z M. 
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[0043] After this, in step 105, the predictor is projected with the aid of a corrector method, which 
is, for example, a numerical Newton method, onto the manifold of the pareto-optimal points. In 
this way a new valid operating point is determined on the manifold of the pareto-optimal points. 

[0044] The steps 103,104 and 105 are repeated iteratively, with the operating point determined 
in the preceding iteration step being used as the starting point for calculating a new valid 
operating point. 

[0045] In step 106 a check is made to determine whether an abort criterion has been met, in 
other words whether, for example, a predetermined number of iterations have been performed 
or a predetermined time limit has been reached. If this is not the case, a return is made to step 
103 and the next iteration is performed. This is continued until the abort criterion is met. 

[0046] Once the abort criterion has been met in step 106, in a next step 107 the set of 
determined pareto-optimal points is restricted to those points in which the weighting factors ctj 

are positive. 

[0047] From these pareto-optimal points, in a final step 108, the user selects an efficient 
operating point of the technical system appropriate to his/her requirements and the technical 
system is designed using this efficient operating point. 

[0048] RfweFig. 2 shows a two-dimensional graphical representation of the predictor-corrector 
method , which may be used in the i nv e nt i on design method . 

[0049] In FigweFig. 2, z' designates a pareto-optimal point on the submanifold M, with this point 
having been obtained in the i-th iteration step of the method. In order to determine a new 
pareto-optimal point, the tangential plane T z jM to the submanifold M is first determined in the 

point z { . The tangential plane is indicated by m e an s of by dashed lines in ftgweFig. 2. In the 

next step 104, a predictor point y i+1 is then determined using normally distributed random 
numbers in the tangential plane T z jM. In the following corrector step 105, which can be, for 

example, a Newton method, the new pareto-optimal point z' +1 is determined. The method is 

then continued, with the pareto-optimal point z' +1 being used as the starting point for new 
predictor step. 
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[0050] FigweFig. 3 relates to a variation of the method accord i ng to th e i nv e nt i on , wherein if 
predictors with negative weighting factors q occur, a reflection is performed in order to 

determine a new predictor with positive otj. RaweFig. 3 shows this projection step being 

performed in a three-dimensional representation. 

[0051] RaweFig. 3 depicts a case in which, starting from a pareto-optimal point z, a predictor 
Yneg is determined which has a negative ctj. This is illustrated graphically in that the section 
between the point z and the point y penetrates the tangential plane T Z M in the point S. In this 
case the point S in turn lies on a subplane of the tangential plane T Z M, for the points of which ' 
the coordinate ctj has the value zero. In order to perform the reflection, the point of intersection 
S is determined first. This can be done using a projection operator which projects the aj 

component from a parameter representation of the straight line running through the points z and 
y. After the point S has been determined, the vector xneg running between S and y can now be 
determined. This vector is then dissected into the tangential component t to the subplane and 
into a normal component n. Thus, t = x neg -n applies to the tangential component. The reflection 

step is then performed, with the new reflected vector x neu having the same tangential 
component t as the old vector x neg and the normal component corresponding to the normal 
component n of the old vector x neg with the sign reversed. The new vector is therefore x neu = t- 
n = (x neg -n)-n = X neg - 2n. There thus results a new predictor y neu which was reflected at the 
tangential plane T Z M. y neu = S + X neu applies to the new predictor point y neu . The above 
described reflection method increases the numerical efficiency of the method since the 
generation of points with negative weighting factors aj is avoided and consequently the 

technical computing resources are used more efficiently. 

[0052] Fig. 4 shows a processor unit PRZE for performing the method accord i ng to th e 
i nv e nt i on . The processor unit PRZE comprises a processor CPU, a memory MEM and an 
input/output interface IOS which is used in a different way via an interface IFC: An output is 
made visible on a monitor MON via a graphical interface and/or output on a printer PRT An 
input is made via a mouse MAS or a keyboard TAST. The processor unit PRZE also has a data 
bus BUS which provides the connection from a memory MEM, the processor CPU and the 
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input/output interface IOS. Additional components such as, for example, additional memory, 
data storage (hard disk) or scanner can also be connected to the data bus BUS. 

[0053] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention covered by the claims which may 
include the phrase "at least one of A, B and C" or a similar phrase as an alternative expression 
that means one or more of A, B and C may be used, contrary to the holding in Superauide v. 
DIRECTV. 69 USPQ2d 1865 (Fed. Cir. 2004V 
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ABSTRACT 

Th e i nv e nt i on r el at e s to| n a method for designing a technical system . Accord i ng to sa i d 
m e thod , a technical system is mod e l l ed modeled by m e an 6 of bv a pre-determined quantity of 
target functions depending on parameters, each individual target function being weighted with a 
weighting factor. The i nv e nt i v e method solves a system of equations comprising the parameters 
and the weighting factors as variables in a variable space, solutions of the system of equations 
forming working points of a solution space in the variable space. SateThe working points are 
determined by m e an 6 of by a predictor-corrector method, according to which a predictor 
produced by m e ans of by a stochastic variable is determined in the variable space, from a first 
working point, and a second working point is then determined in a correcting step. SaidThe 
determined working points are used to design the technical system. The inv e nt i v e method can 
be used to redesign, modify or adapt an already existing technical system. 



